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Perfect fluid spheres, either Newtonian and relativistic, are the first step in developing realistic 
stellar models (or models for fluid planets). Despite the importance of these models, explicit and 
fully general solutions of the perfect fluid constraint in general relativity have only very recently 
been developed. In this Brief Report we present a variant of Lake's algorithm wherein: (1) we 
re-cast the algorithm in terms of variables with a clear physical meaning — the average density and 
the locally measured acceleration due to gravity, (2) we present explicit and fully general formulae 
for the mass profile and pressure profile, and (3) we present an explicit closed-form expression for 
the central pressure. Furthermore we can then use the formalism to easily understand the pattern of 
inter-relationships among many of the previously known exact solutions, and generate several new 
04 , exact solutions. 
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I. INTRODUCTION 



CN 
> _ 

Perfect fluid spheres, both Newtonian and relativistic 0, , have attracted and continue to attract considerable 
attention as the first step in developing realistic stellar models (or models for fluid planets). Whereas some steps 
toward finding all possible solutions to the perfect fluid constraint in the absence of a specific equation of state were 
presented in work early of Wyman and Ho j man et <^ EL 01 > explicit and fully general solutions of the perfect fluid 
' constraint have only very recently been developed 0, Q- In this article we present a variant of Lake's algorithm 
. using curvature coordinates wherein: 



We re-cast the algorithm in terms of variables with a clear physical meaning — the average density and the "grav- 
ity profile" , a quantity closely related to both the gravitational redshift and the locally measured acceleration 



^jy due to gravity. 



• We minimize the number of differentiations and integrations by several judicious applications of integration by 
parts. 



• We present explicit, compact, and fully general formulae for the mass profile and pressure profile of an arbitrary 



fluid sphere. 

• We present an explicit, compact, and general formula for the central pressure of an arbitrary fluid sphere. 

• We compare and contrast the relativistic formulae wc obtain with the much simpler Newtonian situation. 

We emphasise that one of the virtues of this type of approach is that one is not fixed a priori to dealing with a 
pre-specified equation of state [g — in many interesting physical situations the equation of state is either uncertain 
or, because the fluid in question might be inhomogeneous, it may not even make sense to assign a single equation of 
state to the entire fluid sphere. 

To further illustrate the formalism we show how it may be used as the basis for a partial classification scheme 
— there is a free parameter in the algorithm that can take simple solutions into more complicated ones. Once this 
is appreciated it becomes easy to see (simply by parameter counting) that certain simple solutions must have one- 
parameter generalizations. Conversely, this observation explains why so many of the earliest discovered exact solutions 
have one-parameter extensions. 
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II. FRAMEWORK 



To set the stage, consider a static spherically symmetric geometry. It is a standard result that without loss of 
generality we can choose coordinates to write the metric in the form 



ds 2 



exp 



g(f) df 



dt 2 



dr 2 



1 — 2m(r)/r 

Here g{r) "gravity profile" . It is related to the gravitational redshift by 



r 2 [d9 2 + sin 2 d<f> 2 ] 



1 + z = exp 



g(f)df 



and is related to the locally measured acceleration due to gravity by 



2m(r) 

i g(r)- 



(1) 



(2) 



(3) 



Our convention is that g(r) is positive for a downward acceleration. The function m(r) is the quasi-local mass. 
In the vacuum region beyond the surface (if any) of the star-like object, the Schwarzschild solution yields g(r) = 
(M/r 2 )/(l — 2M/r) and m(r) = M. We find it more convenient to write the metric in the form 



ds 2 



exp 



](f) df 



dt 2 



dr 2 



1 - 2/i(r) r 2 



r 2 [d<9 2 + sin 2 9 d</> 2 ] 



(4) 



where p(r) — ^-p{r) is proportional to the average density inside radius r. In terms of these variables, the Einstein 
equations are 

8np = G u = 2m'{r)/r 2 = 2[r p!{r) + 3/z(r)]; 



_ Jff(r) [l-2//(r)y 



87175 = Gff = 2 



8np = G m = -r[l + rg(r) 



dn(r) 



-MO ; 



2^[l + r 3 (r)] 2 +r 2 ^ 



dr [ dr 

The first of these equations integrates to 

1 r 

A*( r ) = ~5" / 47rp(f) f 2 df 
r Jo 

which justifies the choice of notation m(r) = p(r) r 3 . 



dg(r) g(r) , 
J ^ + — +g{r) 2 

dr r 



(5) 
(6) 

(7) 
(8) 



III. GENERAL SOLUTION AND GENERATING FUNCTION 

By demanding the isotropy condition Gff. = G^ and algebraically solving for dg/dr we obtain 



dg 
dr 



'„3 



= -9 



1 + n'r 



9 



r /i 



r(l - 2(j, r 2 ) * 1 - 2/i r 



(9) 



This is a Riccati equation, for which there is no general analytic solution. If on the other hand we take this same 
equation and rearrange it algebraically to extract dp/dr we find 



d/i _ 2r (g 2 + g') | {g/r)' + g 2 /r 
dr 1 + r g l + r g 



(10) 



But this is a simple first-order linear ODE and hence explicitly solvable. A symbolic manipulation program such as 
Maple, or a slightly tedious hand-computation, easily yields the general solution 



fi{r) = exp 



/ r [g 2 {r)+g'{r)] ^ 
1 + rg(r) 



, -g(r) +rg'(r) +rg(r) 2 
Cl + 1 r 2 [l + rg(r)} 6XP 



-2 / r [g2(r) + < (r)] dr 
1 + rg{r) 



■ (11) 
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This statement is equivalent to the algorithm presented by Lake 0: Given a prescribed gravity profile g(r) (the 
"generating function"), and the knowledge that we are dealing with a perfect fluid, the mass profile m(r) = fi(r) r 3 
is deduced in closed form. The algorithm (|lf (I is also equivalent to that presented by Rahman and Visser 0, after a 
change of coordinates (from isotropic to curvature coordinates) and a change of variables. A particularly nice feature 
of the present version of the algorithm is that the generating function g(r) has a clear physical interpretation in terms 
of the gravitational field. Now the above is by no means the most useful form in which fj,(r) can be presented. An 
integration by parts permits us to simplify the appearance of the integrating factor 



exp 



r [9 2 {r)+g'(r)} 



dr 



1 + rg{r) 

It is now extremely useful to introduce the notation 



[1 + r g(r)] 2 exp 



9(r) 



1 - r g{r) 
1 + r g(r) 



dr 



0(r) 



1 -r g{r) 

din — — dr; 

1 + r g[r) 



and 



■d{r 1 \r 2 ) 



1 - r g(r) 

9 \r) — r^r dr. 

1 + r g(r) 



(12) 



(13) 



We warn the reader that we cannot generally assume r g(r) < 1, and the consequently $ may become negative. For 
instance, in the physically reasonable regime m(r)/r > 1/3 [that is, /i(r) r 2 > 1/3] and p > 0, it can be shown from 
the Gff Einstein equation that r g(r) > 1. All that we can safely say in general is that as long as local gravity points 
down wc must have 



— J g(r) dr < i) < J g(r) dr 



With this notation 



fi(r) 



exp[+2??(r)] 



C 2 



[1 + r g(r)] 2 
A second integration by parts now yields 
exp[+2tf(r)] 



[1+r <?(r)] 



~$( r ) + r 9'( r ) + Tg{rY 



exp[-2i?(r)] dr 



[l + r g(r)f 



C 3 



, ^2 exp[-2^(? 
r 9{r)} -o 



9(r) 



exp[-2i?(r)] dr 



(14) 



(15) 



(16) 



In this version of the result we have eliminated all the derivatives of g(r). A third integration by parts, using 



1/r 3 



i(l/r 2 )' then leads to 



(j,(r) 



g(r) [l + \rg{r)\ exp[+2^(r)] 



r [1 + r g(r)] 2 



■ g{r)f 



C 4 + 2 



9{rf 



r[l +r g(r)} 



exp[-2i?(r)] dr 



(17) 



This final version, as we shall soon see, has nice behaviour at the origin. Again, we emphasise that this is the explicit, 
and most general, solution to the perfect fluid constraint for arbitrary generating function g(r). All perfect fluid 
spheres, no matter how derived, must satisfy this equation. 

The pressure can now be determined using the Gff Einstein equation so that 



p(r) = 



1 



87r[l + r g(r)} 2 



-g(r) 2 - 2[1 + 2r g(r)} exp[+2tf (r)] <C 4 + 2 



9{rf 



r[l + r g(r)\ 



exp[-2i?(r)] dr 



(18) 



This now provides for us the explicit and fully general solution to the mass profile and pressure profile, given only 
the gravity profile and the information that we are dealing with a static spherically symmetric perfect fluid. For a 
consistency check, we can compare these formulae to the much simpler result for Newtonian stars: 



/i(r) 



flO) 

r 



and 



P(r) 



1 



-9{r? + c 5 



9{r? 



dr 



(19) 



To complete the analysis we should now impose boundary conditions. There are three natural locations to work with: 
(1) the center of the fluid body, (2) the surface of the fluid body [assuming it has a well defined surface], and (3) 
spatial infinity. Perhaps surprisingly, the simplest results are obtained if we normalize at spatial infinity. 
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IV. BOUNDARY CONDITIONS AT SPATIAL INFINITY 



We will now adopt the very mild condition that the total mass of the fluid sphere is finite, so that /i(r) — > as one 
approaches spatial infinity. We also assume p(r) — > at spatial infinity. Then (from the Gff equation) we deduce 
g(r) — > at spatial infinity. Physically this means that the present discussion is capable of handling situations with 
a tenuous atmosphere extending all the way to infinity, and that the special case where the fluid body has a sharp 
surface with p(r > R) = and m(r > R) = M is automatically included. Then fixing boundary conditions at spatial 
infinity, the mass profile If 7(1 is given by 



n _ffW [1 + |r 9(r)] _ exp[-2tf(r; oo)] 



9{f) 



,-'\2 



r [l + rg(r)] 2 [l + rg(r)] 2 

We can simplify this slightly to yield 

, * = [l + hrg(r)} _ 2 

r [l + rg(r)} 2 [1 + r g(r)} 2 J r f[l + fg(f) 

The pressure profile determined from If 8|) is then 



f[l + r g(r)] 
5« 2 



exp[+2i?(f ; oo)] df . 



p(r) 



f 



8tt[1 + r g(r)} 2 
For the central pressure, p c = p(0), we find 



-g(r) 2 +4[l + 2r g(r)} 



f[l + f g(r)] 



exp[— 2i?(r; f)] df. 



exp[— 2i?(r; f)] df 



g(r) 2 



exp[— 1?(0; r)] dr. 



2tt 7 r[l + r g(r)} 

Compare with the equivalent statement for a Newtonian fluid body in which the pressure profile is 



p{r) 



1 

8^ 



and the central pressure is 



~9{rf 



2tt 



9{r) 



g{rf 



dr. 



df 



(20) 



(21) 



(22) 



(23) 



(24) 



(25) 



V. BOUNDARY CONDITIONS AT THE CENTER OF THE FLUID SPHERE 



If we apply boundary conditions at the center of the sphere then, using 

m(r) + 4irp(r)r 3 



g(r) 



r 2 [l — 2m(r)/r] 



and the assumed finiteness of p c and p c , implies 



47T 



g(r) = —(p c + Zpc)r + 0{r z ). 
The mass and pressure profiles are given by 



r [l + rg(r)} 2 [l + rg(r)] 2 



9(f) 2 



o r[l + rg(r) 



exp[-2i?(0; f)] df 



and 
p(r) 



I 



8?r[l + r g(r)} 2 
to be compared with the Newtonian result 



g{rf + [1 + 2r g(r)} exp[+2tf (0; r; )] <^ 8tt p c - 4 



g(ff 



(26) 



(27) 



(28) 



o f[l + fg(f)] 



exp[-2z?(0; f)] df 



P(r) =Pc~ 



8tt 



g(rY + 4 



(29) 
(30) 
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VI. BOUNDARY CONDITIONS AT THE SURFACE OF THE FLUID SPHERE 



If the fluid sphere has a sharp boundary, (say at radius R, with the density and pressure identically zero outside 
this radius), then it can be useful to normalize at this surface. For the pressure profile we find [in terms of the "surface 
gravity" g s = {M/R 2 )/^l - 2M/R] that 



p(r) 



1 



87r[l + r g{r)] 2 



-g(r) 2 + [1 + 2r g(r)] exp[-2tf(r; R)} 



9(r) 2 



l + 2Rg s J r f[l + fg(f)] 



exp[+2tf(f; R)} df 



There is a similar but uninteresting expression for /i(r). The central pressure is now 

if g gexp[-2tf(0;fl)] 4 [* 9ir? A 
Pc -8n\ l + 2Rg s + A J r[l + r g(r)] 6Xp[ ^ (U ' ^ ' 



(31) 



(32) 



where we have now reduced the range of integration from (0, oo) to (0, R) at the price of introducing an extra term 
depending explicitly on total mass and radius of the fluid sphere. This can be compared with the equivalent Newtonian 
results 



p(r) = 



1 

8^ 



,2 „f„\2 i a I v\ i ^ 



gi - g(ry + 4 



and 



with g s = 



M 
R 2 ' 



(33) 



(34) 



VII. SOLUTION GENERALIZATION TECHNIQUE 



One nice feature of the present analysis is that it allows one to turn simple exact solutions into more complicated 
ones. While the general algorithm presented above always provides an exact solution, it may not be an "elementary" 
solution in the sense that the integrations might not be do-able in terms of either elementary or special functions. In 
such a situation, a simplified algorithm is sometimes useful. 

Suppose that one has found, by some unspecified means, a specific exact solution for a perfect fluid sphere. Let 
that exact solution be given in terms of m(r) [or equivalently /i(r)] and g(r). Then for any arbitrary constant fc, 



/i(r) — > fi(r) + k exp 



-2 



r [9 2 (r)+g'{r)] 
1 + rg(r) 



dr 



(35) 



is also an exact solution for a perfect fluid sphere [with the same g(r)]. This construction may sometimes "fail" in the 
sense that the integral is either too trivial [returning you to the seed solution you started with] , or too complicated to 
perform in terms of elementary or special functions. However in very many cases this simple construction is sufficient 
to understand why certain broad classes of exact solution exist. 

Let us start by rescaling the time variable to remove any redundancies in the number of free parameters, n, 
appearing in g u ■ If the number of free parameters appearing in g rr is not at least n + 1 then the seed solution you 
have must have a generalization. For instance the Minkowski solution [a particularly simple fluid sphere with zero 
pressure and density] has exactly zero parameters appearing in both g tt and g rr , and so must have a one-parameter 
generalization. In this case, performing the integration leads to the Einstein static universe. Similarly, the exterior 
Schwarzschild solution [another particularly simple fluid sphere with zero pressure and density] has exactly one free 
parameter [the mass] appearing in both g u and g rr , and so must have a one-parameter generalization. In this case, 
performing the integration leads to what is called the Kuch68 II solution in the Delgaty-Lake classification A 
slightly more complex example, using anti-de Sitter space as a seed, leads to the Tolman IV solution. A number of 
additional examples of this phenomena are collected in table I. 

Of course, sometimes explicit exact solutions were first discovered in their general form, in which case this algorithm 
provides no extra information. (This comment applies, for instance to the Wyman lib geometry.) Conversely, 
sometimes the integral is too complicated to provide a closed-form solution — the generalization my be exact but 
too complex to write down explicitly. (As for instance when you use the Schwarzschild-de Sitter [Kottler] geometry 
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Seed 


Generalization 


Minkowski 


Einstein static 


Schwarzschild exterior 


Kuch68 II 


anti-de Sitter 


Tolman IV 


Tolman V 


Kuch2 I 


Tolman VI 


Wyman Ha 


Kuchl lb 


appears new 


M-W III 


appears new 


K-0 III 


appears new 



TABLE I: Table I: Seed solutions and their generalizations. 



as seed. Similarly, by parameter counting Tolman VII and Tolman VIII must have one-parameter extensions, but it 
seems impossible to write then down in closed form.) 

Finally, we point out that there are some cases where this formalism does lead to apparently new solutions. (We 
again follow the Delgaty-Lake classification Q-) For instance, the Kuchl lb solution 



dr 2 

ds 2 = -(A r + B rlnrfdt 2 + — + r 2 dft 2 (36) 



generalizes to 



d.s 2 = (A r + B rlnrfdt 2 + - f \ + ™^ r) + ^ 2 dr» + r 2 dtt 2 (37) 
v ' (2A + 2Bln(r) + B) - kr 2 y ' 



which appears to be new. Similarly, the M-W III solution, which can be cast into the form, 

ds 2 = - (r - -) dt 2 + - 7dr ' +r 2 dn 2 (38) 
\ a J 4(1 - r 2 /a 2 ) 

generalizes to 

*r=\r--) dr + — dr 

ith 



ds'=[r--) dt 2 + +r 2 dn 2 (39) 

1 a J 1 — 2m{r)/r 



. . 4r 2 + 3a 2 , (r - a)r 10 / 3 



which also appears to be new. Also, the K-0 III solution can be cast into the form 



ds 2 = - ( 1 + ^ ) dt 2 + dr 2 + r 2 dfl 2 (41) 



which is spatially flat. It generalizes to 



a 2 



,2 \ 2 ,2 



dg2 = -{ 1 + ^) ^ U-kr 2 ^^ 2 )- 2 ^ ^ 2 ^ (42) 
which is contained within the new class of exact solutions briefly described by Lake . 



VIII. DISCUSSION 



As emphasised in the article by Rahman and Visser |6j, and reiterated by Lake jjj, while this type of algorithm 
guarantees a perfect fluid body it does not necessarily guarantee a "physically reasonable" perfect fluid body. One 
physically reasonable constraint that is easy to enforce in the current formulation is g > 0; locally measured gravity 
should always attract towards the center of the body. A second physically reasonable constraint which is automatically 
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satisfied is that the central pressure is positive. It is considerably more difficult to enforce m(r) > 0, p(r) > 0, and 
p(r) > 0. Checking these physically motivated constraints amounts to mathematically investigating a set of integral 
inequalities, and seems to require a case by case investigation depending on the assumed gravity profile g(r). One 
should not however lose track of the significance of what has been accomplished: 

• We have derived the exact and fully general solution to the pressure isotropy condition in terms of variables 
that have a direct physical meaning, the gravity profile g{r) and mass profile m(r). 

• We have also derived an exact and fully general formula for the pressure profile p(r) of a perfect fluid sphere 
that depends only on the gravity profile g(r). 

• In particular we have an exact and fully general expression for the central pressure of a fluid sphere, again 
determined directly in terms of the gravity profile g(r). 

• The algorithm provides a natural framework for understanding the reason for the existence of certain broad 
classes of exact solution, and in some cases leads to new exact solutions. 

Because this algorithmic approach works directly in terms of physically meaningful quantities, with a physically 
meaningful "generating function" in the form of the gravity profile g{r), the interpretation of the results is somewhat 
clearer than in the algorithms presented in the Rahman- Visser [(| and Lake Q articles. We expect that this version 
of the algorithm for generating perfect fluid spheres will lead to additional useful "exact solutions" . In particular, the 
new class of exact solutions briefly described in has a very natural representation in terms of this algorithm. 
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